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Quasar absorption spectra at 21-cm and UV rest-wavelengths are used to estimate the time 
variation of x = a 2 g p /i, where a is the fine structure constant, g p the proton g factor, and m e /m p = 
H the electron/proton mass ratio. Over a redshift range 0.24 < z a b s ^ 2.04, (Ai/i)2 llted = 
(1.17 ± 1.01) x 10 -5 . A linear fit gives x/x = (-1.43 ± 1.27) x 10~ 15 yr~\ Two previous results on 
varying a yield the strong limits Afi/fi = (2.31 ± 1.03) x 10" 5 and Afi/fi = (1.29 ± 1.01) x 10" 5 . 
Our sample, 8x larger than any previous, provides the first direct estimate of the intrinsic 21-cm 
and UV velocity differences ~ 6 km s _1 . 

PACS numbers: 98.80.Es, 06.20.Jr, 95.30.Dr, 95.30.Sf 
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The existence of extra spatial dimensions, often in- 
voked by superunification theories, may be inferred by 
the detection of spatial or temporal variations in the val- 
ues of coupling constants (See 0] for a review). Spec- 
troscopy of gas clouds which intersect the lines of sight 
to distant quasars is a unique tool, probing the values 
of these constants over a large fraction of the age of the 
universe. The highly sensitive many- mult iplet method, 
developed by 0, 01 > nas been applied to rest-frame ul- 
traviolet (UV) atomic quasar absorption lines to provide 
constraints on the possible variation of the fine struc- 
ture constant, a = e 2 /(hc) 0, 0, H H 0- Molecular 
hydrogen absorption lines have provided constraints on 
the variability of the electron-to-proton mass ratio, fi 

Principle. — The above results have been obtained 
by use either of heavy element transitions, which ab- 
sorb in the rest-frame UV, or of molecular hydrogen 
transitions. Another approach is to use the parameter 
x = a 2 g p m e /m p when, as well as rest-frame UV, rest- 
frame 21-cm absorption, due to cold neutral hydrogen, is 
also detected. Rest-frame UV absorption is observed red- 
shifted in the optical region, and rest-frame 21-cm is ob- 
served redshifted at longer radio wavelengths. The ratio 
of frequencies a>2i/^uv oc x. 21-cm absorption occurs in 
a few damped Lyman a (DLA) systems which also show 
heavy-element absorption in the UV. A detailed list can 
be found in If both UV and 21-cm absorption oc- 
cur at the same physical location, the relative change of 
the value of x between redshifts z and is related to the 
observed absorption redshifts for rest-frame 21-cm and 
UV, z 2 i and z uv , according to Ax/x = = Z ^~ Z21 

I 1' Z1 U V ' ° ' X 1 + 221 

|14| . We obtained values for z\jy and Z21 by using the 
strongest absorption components in an absorption sys- 
tem. This approach is discussed in detail later. 

However, as there are only 17 DLAs where both 21- 
cm and UV absorption have been detected, there are 
few results based on this method [3 [H, [3, Out 



of these, only 0] use high-resolution optical data from 
the Keck telescope's HIRES spectrograph, but they pro- 
vide an estimate of Ax/x at a single redshift from a sin- 
gle absorption system. We applied this method to eight 
absorption systems in eight quasar spectra (one system 
per spectrum), covering the absorption redshift range 
~ 0.24 to ~ 2.04. We used all available 21-cm absorption 
data in conjunction with the highest-resolution UV data 
available. Thus the results presented here are based on 
the largest dataset of the highest quality to which this 
method has been applied to date. 

Data analysis. — Details of the 21-cm and UV data 
used are given in Table [I] All redshifts are in the helio- 
centric frame. For the strongest component in each 21- 
cm absorption complex, the dispersion coordinate at the 
pixel of minimum intensity, MHz or km s _1 , was mea- 
sured, from which Z21 was obtained. We searched the 
optical data for heavy element absorption features close 
to the redshifts where there is 21-cm absorption. A num- 
ber of UV absorption features were thus identified, some 
due to neutral species and most due to singly ionized 
species. For all UV spectra we determined the value of 
the dispersion coordinate, A or km s _1 , for the strongest 
component at the pixel of minimum intensity. We then 
determined absorption redshifts for each detected neu- 
tral or singly ionized absorption species that was not 
saturated. A zuv value was determined individually for 
each transition of a single species, e.g. independently for 
Znn 2026.14 and Znn 2062.66. In all, there were 30 dis- 
tinct UV species identifications (see Table H] column 6). 
Detailed velocity plots showing all 21-cm and UV absorp- 
tion components used can be found in 

[3 and [13. 

Estimating Ax/x. — For each absorption system we 
calculated (zuv), the average of all UV absorption red- 
shifts for single species (column 5 in Table QJ. Using 
this and our measured Z21 (column 3 in Table we ap- 
plied the relation between Ax/x, zjjy and Z21 to obtain 
(Ax/x)ab s . This is the average value for Ax/x for each 
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TABLE I: Data used in this work. There is one 21-cm/UV absorption system in each quasar spectrum. Column 1 is the 
quasar name and Column 2 its emission redshift. Column 3 gives the 21-cm absorption redshift (and error from Column 4 
references) for the strongest component. We determined this after digitizing 21-cm absorption plots (references in column 4). 
For Q0458— 0203 the original data were used with the error taken from |2C| . Column 5 gives the mean absorption redshift (and 
standard deviation on the observed mean) for the strongest UV component. Column 6 gives the UV heavy element species 
observed in the optical (with number of transitions in parentheses, if more than one). Column 7 gives the source for the UV 
data. Data for eight quasars were obtained from the European Southern Observatory's (ESO) archive and were originally 
observed with the UVES spectrograph on the Very Large Telescope (VLT), in which case the ESO program ID is given in 
column 7 and the principal investigators of the program are given in footnotes. For Q1331+170 we also used Sill 1808.01A 
Keck/HIRES data provided by A. Wolfe. For quasar Q0235+164 we digitized an absorption plot from the literature for a single 
heavy element species. 



Quasar 



Z21 



21-cm data 



(guv) 



UV data 



Q0952+179 
Q1127-145 
Q1229-021 
Q0235+164 
Q0827+243 
Q1331+170 



1.478 0.237803(20) 

1.187 0.312656(50) 

1.038 0.394971(4) 

0.940 0.523874(100) 

0.941 0.524757(50) 

2.097 1.776427(20) 



Q1157+014 1.986 1.943641(10) 
Q0458-0203 2.286 2.039395(80) 



[21] 0.237818(6) Mgi, Caii(2) 

[22] 0.312648(6) Call(2), Mnn(3) 

[23] 0.395019(40) Caii(2), Mnii(3), Tin 

[24] 0.523829(6) Mgi 

[21J 0.524761(6) Caii(2), Fen 

[22] 1.776355(5) Mgi, Al 11, Sin, S 11, 

Cl(3), Cl*, Crn(2), Mnn(2), 
Fen (4), Nin(6), Znn 
Sin c 

[1JJ 1.943738(3) Mgi, Mgn(2), Sin, 
Nin(6) 

[22] 2.039553(4) Znii(2), Nin(6), Mnii(3), 
Cr n(3) 



69.A-0371(A) a 

67. A-0567(A) 6 , 69.A-0371(A) a 

68. A-0170(A) C 
[25] 

68.A-0170(A) C , 69.A-0371(A) a 
67.A-0022(A) d , 68.A-0170(A) C 
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65.O-0063(B) / , 67.A-0078(A) / , 
68.A-0461(A) 9 

66.A-0624(Ay, 68.A-0600(A) / , 
072.A-0346(A) / , 074.B-0358(A)'' 
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absorption system. We plot these results in Fig. ^ 

Taking the average of all {Ax/x) a hs values, we ob- 
tained (Ai/x) total = (0.91 ± 1.04) x 10~ 5 ('result 1') 
over an absorption redshift range 0.24 < z a b s < 2.04 
and a fractional lookback time range 0.20 < tflb ^ 0.76, 
where we have used a Hubble parameter Hq — 73 km 
s _1 Mpc -1 , a total matter density Om = 0.27 and a cos- 
mological constant — 0.73. The error quoted is the 
standard deviation on the mean. 

We also performed the above calculations by taking 
into account statistical errors on Z21 and (ztjv) P er aD - 
sorption system (Table[IJ. For Z21 we used errors from the 
Table references |32|. For (ztjv) we used the standard 
deviation on (ztjv) for each absorption system [33}. We 
obtained {Ax/x)™"^ tcd = (2.18 ±0.97) x 10" 5 ('result 
2'). This has a \ 2 P er degree of freedom, v, x 2 l v = 8. 
We thus increased the individual errors on (Ax/x) a b s 
by an additional error, s, to (& 2 Ax / x \ b + s 2 ) ' 5 until, 
at s = 1.90 x 10 -5 , x 2 l v = 1. This corresponds to 
{Ax/x}™^ tcd = (1.17 ± 1.01) x lO" 5 ('result 3'). 

For all (ztjv) values per quasar absorber we calcu- 
lated the average fractional lookback time per absorber, 
(^flb)abs- We then performed an iterative linear least 



squares fit to (Ax/x) a ^ s = A(£flb)abs, where the ad- 
ditional error, s, was determined at each iteration to 
force x 2 l v — 1 around the fit, obtaining A — (1.90 ± 
1.69) x 10- 5 [13- lt follows that the best fit rate of 
change of x as a function of time 4z(Ax/x) = x/xq = 
(-1.43 ± 1.27) x 10- 15 yr -1 - 

Assumptions.- In this work we are making two as- 
sumptions. The validity of our result does not depend on 
the validity of these assumptions. On the contrary, we 
are essentially using our result to test these assumptions. 
Moreover, if these assumptions are incorrect, they only 
contribute to any observed scatter in (Acc/a:) a bs- There- 
fore, this possibility has already been taken into account 
in results 1 and 3. 

1. Strongest components: Both the 21-cm and UV pro- 
files exhibit complex velocity structure, i.e. have multi- 
ple absorption components at slightly different redshifts. 
That being the case, how does one compare the red- 
shifts among different transitions? For neutral and singly 
ionized UV species, velocity structure is essentially the 
same and corresponding components can easily be iden- 
tified. This is not the case if one compares 21-cm and UV 
velocity structure, although we have not systematically 
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Ax/x = (1.17±1.01) x 10" 
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FIG. 1: Ax/x results for the eight abosrption systems in 
our quasar sample. Initial (increased) error bars have shorter 
(longer) terminators. Quasar names are given truncated to 
four digits. Each point represents (Ax/x) a b s obtained from 
Z21 and (zuv), for all heavy element species in a quasar spec- 
trum, versus average (zuv) for that spectrum. The solid hor- 
izontal line is result 3. The dashed lines show the ±1<t range. 



searched for corresponding velocity patterns in the 21-cm 
and UV profiles. One simple option is to measure mean 
absorption centroids over all components for each absorp- 
tion system, which was essentially the approach taken by 
|l6| . Alternatively one can simply determine Z21 for the 
strongest component in the 21-cm profile and ztjv for the 
strongest component in each of the UV transitions, and 
use these. In our sample this is the only well defined 
quantity for both the 21-cm and the UV profiles Addi- 
tionally, this procedure, unlike mean centroiding, should 
have the advantage of being less sensitive to measurement 
errors caused by broad velocity structure in absorption 
lines. 

2. Neutral and singly ionized species: We have used 
both neutral and singly ionized species. One might ex- 
pect only the UV neutral gases to be spatially coincident 
with the 21-cm absorbing gas (unless the 21-cm gas is 
primordial or has very low heavy clement abundances). 
However, it was clear in our sample that velocity struc- 
ture of neutral UV species was followed closely by singly 
ionized species as well. This suggests that ionization frac- 
tion and abundance may not change significantly along 
a complex. As shown by [13], this is a general obser- 
vation for DLAs. Further, the strongest component in 
singly ionized species corresponded well to the strongest 
component in neutral species. These observations justify 



using singly ionized species as well. As neutral species 
are very rare, one obtains the substantial advantage of a 
much larger sample. This allowed us to investigate sys- 
tematics for individual absorbers for the first time (see 
below) . 

Errors. — The redshift errors introduced in digitiz- 
ing the data are < 4 x 1CP 6 , (typically ~ 10~ 6 ), which 
is about an order of magnitude smaller than the ob- 
served intrinsic scatter on Ax/x. Further, for result 1, 
the error quoted is a standard deviation on the mean, 
directly reflecting the scatter in (Ax/x}nbs- The error 
in result 2 was obtained after introducing errors for Z21 
and (ztjv) (Table QJ. Although, formally, this gives a 
smaller error, when the errors in (Ax/a;) a bs are increased 
so that x 2 / 1 ' = U results 1 and 3 for (Ax/x)™°^ tcd 
and its error are very close. This shows that the scatter 
in (Ax/x) a bs dominates any errors in individual redshift 
measurements. For a single absorption complex this is 
a systematic error. However, for many complexes this 
error has random sign and magnitude. This error can 
be accounted for directly if treated as a statistical error 
when calculating the mean value of x for many absorp- 
tion complexes, as we do for results 1 and 3. 

Discussion of results. — In Fig.^there is considerable 
scatter in the values for (Ax/x)^. For all spectra all op- 
tically observed species tend to group together in their 
Ax/x values on one side of zero. This is because there is 
significant offset between the single Z21 value and all ztjv 
values in a system. This suggests that there is spatial 
offset between the 21-cm and UV absorbing gases, but, 
as expected, this is random for different absorbers. The 
error of result 1 directly reflects this scatter. A straight- 
forward statistical calculation misses this effect (small 
error in result 2). Result 3 is consistent with result 1 and 
the value of s provides an estimate of the line of sight 
velocity difference Awi os ~ cs = 6 km s . 

One possible physical explanation for the observed off- 
set may be a large angular size for the emitting 21-cm 
quasar source, as seen by the absorber. A 21-cm sightlinc 
can then intersect a cold, neutral hydrogen cloud with 
little or no heavy elements, whilst a UV/optical sight- 
line can intersect another cloud with heavy elements at 
quite a different velocity. A large angular size is due to 
the combined effects of proximity of the absorber to the 
quasar (small z cm — (z)) and intrinsic size of the radio 
emitting region. Note also the good agreement with the 
velocity difference (~ 10 km s _1 ) between 21-cm and mm 
absorption lines likely due to small scale motion of the 
interstellar medium [28| . 

Comparison with previous results. — used neu- 

tral carbon C 1, Ci*, lines in the Keck/HIRES absorp- 
tion spectrum of Q1331+170, to obtain Ax/x = (0.70 ± 
0.55) x 10~ 5 . For the same object we obtain Ax/x — 
(-2.59 ± 0.74 stat ± 1.90 syst ) x 10~ 5 . The central value in 
[16j differs from ours because these authors used a ztjv 
value which is a weighted mean for the observed compo- 
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nents 29], thus obtaining a value very close to z 2 \ for 
this absorption system. In our VLT/UVES spectrum of 
Q1331+170, we detected 23 distinct UV heavy element 
transitions whose strongest component was well defined 
within a few km s . Further, our use of eight objects 
allows us to quantify the systematics due to the sightlinc 
issues explained above. The error in 16] is exclusively 
statistical, as it is based on a single absorption system. 
Although at face value this error is lower than ours, it in- 
evitably contains no information on systematics, which, 
as our error estimate shows, dominate. 

Robustness. — We stress that a non-zero Ax/x value 
is not corroborated by the sample as a whole, for which 
the result is robust. For our quoted result we have not 
used zyv obtained from Ca n whose ionization poten- 
tials are least similar to those for all other elements used. 
Even so, the result changes by about 4% if Can is in- 
cluded. If we use Z21 from the literature, rather than 
values determined from our digitized plots, we obtain 
(Ax/x)^ tcd = (1.51±1.04)xlO- 5 H3. If, additionally, 
we do not use results for Q0235+164 or Q1127-145, thus 
excluding any results based on digitized plots, we obtain 
(Ax/x)^ tcd = (1.77 ± 1.12) x 1(T 5 . In all cases Aui OB 
remains at 6 km s _1 pjfij ]. 

Variation of [i. — Measurements of x = a 2 g p m e /m p 
are sensitive to variation of several fundamental con- 
stants. The proton mass m p is proportional to the quan- 
tum chromodynamic (QCD) scale Aqcd- The depen- 
dence of the proton mass on the quar k mass m q is very 
weak (Am p /m p w 0.05Am 9 /m g [30j) and may be ne- 
glected. The dependence of the proton magnetic g p fac- 
tor on the fundamental constants is also quite weak: 

A 5p / 5p w S3" Therefore > the most 

important effects are due to the variations of a and 
fi = m e /m p w 0.2m e /AQCD- From the definition of x 
one obtains ^ A " 1 = Afi/fi = Ax/x - 2Aa/a - Ag p /g p . 
Using the Act /a results of H and 0, we obtain A^/^j, = 
(2.31±1.03)xlCT 5 and Afi/fi= (1.29 ± 1.01) X 10" 5 , re- 
spectively, assuming Ag p /g p ~ 0. Our result contradicts 
Afi/fi = (-2.97 ± 0.74) x 10" 5 pu||37|but is consistent 
with Afi/fi = (0.5 ±3.6) x 10~^ (2er)|ljj. Note that our 
result on Afi/fi is derived using a completely independent 
method compared to 0, IjHJ . 

A more sophisticated analysis, involving fitting Voigt 
profiles, would not significantly affect our results because 
the uncertainties in our estimates of the UV redshifts are 
small compared to the offset between z 2 i and zuv ■ While 
more data will improve our estimate of A«i os , a sample of 
~ 100 21-cm/UV absorbers is required before averaged 
line-of-sight velocity differences and individual redshift 
measurement errors are similar. 

We thank A. Wolfe for data and useful discussions. 
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FIG. 2: Velocity plot for 21-cm and UV absorption towards quasar Q0952+179. The solid vertical line at km s _1 is at 221- 
The dotted vertical line is at (ztjv)- In this and subsequent plots Call is shown for illustration only but has not been used in 
the calculation of the plotted (zuv)- 
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FIG. 3: Velocity plot for 21-cm and UV absorption towards quasar Q1127— 145. The solid vertical line at km s 1 is at 221- 
The dotted vertical line is at (zuv). 
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FIG. 4: Velocity plot for 21-cm and UV absorption towards quasar Q1229— 021. The solid vertical line at km s 1 is at Z2i- 
The dotted vertical line is at (zuv)- 
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FIG. 5: Velocity plot for 21-cm and UV absorption towards quasar Q0235+164. The solid vertical line at km s 1 is at 231. 
The dotted vertical line is at (zuv). 
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FIG. 6: Velocity plot for 21-cm and UV absorption towards quasar Q0827+243. The solid vertical line at km s _1 is at Z2i- 
The dotted vertical line is at (zuv). 
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FIG. 7: Velocity plot for 21-cm and UV absorption towards quasar Q1331+170. The solid vertical line at km s _1 is at 221. 
The dotted vertical line is at (zuv)- 
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FIG. 8: Velocity plot for 21-cm and UV absorption towards quasar Q1157+014. The solid vertical line at km s 1 is at 221- 
The dotted vertical line is at (ztjv). 
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FIG. 9: Velocity plot for 21-cm and UV absorption towards quasar Q0458— 0203. The solid vertical line at km s is at 221 . 
The dotted vertical line is at (zuv). 



